A TCF4-dependent gene regulatory network confers resistance to immunotherapy in melanoma.

To better understand intrinsic resistance to immune checkpoint blockade (ICB), we established a comprehensive view of the cellular architecture of the treatment-naive melanoma ecosystem and studied its evolution under ICB. Using single-cell, spatial multi-omics, we showed that the tumor microenvironment promotes the emergence of a complex melanoma transcriptomic landscape. Melanoma cells harboring a mesenchymal-like (MES) state, a population known to confer resistance to targeted therapy, were significantly enriched in early on-treatment biopsies from non-responders to ICB. TCF4 serves as the hub of this landscape by being a master regulator of the MES signature and a suppressor of the melanocytic and antigen presentation transcriptional programs. Targeting TCF4 genetically or pharmacologically, using a bromodomain inhibitor, increased immunogenicity and sensitivity of MES cells to ICB and targeted therapy. We thereby uncovered a TCF4-dependent regulatory network that orchestrates multiple transcriptional programs and contributes to resistance to both targeted therapy and ICB in melanoma.

A cellular hierarchy in melanoma uncouples growth and metastasis.

Although melanoma is notorious for its high degree of heterogeneity and plasticity1,2, the origin and magnitude of cell-state diversity remains poorly understood. Equally, it is unclear whether growth and metastatic dissemination are supported by overlapping or distinct melanoma subpopulations. Here, by combining mouse genetics, single-cell and spatial transcriptomics, lineage tracing and quantitative modelling, we provide evidence of a hierarchical model of tumour growth that mirrors the cellular and molecular logic underlying the cell-fate specification and differentiation of the embryonic neural crest. We show that tumorigenic competence is associated with a spatially localized perivascular niche, a phenotype acquired through an intercellular communication pathway established by endothelial cells. Consistent with a model in which only a fraction of cells are fated to fuel growth, temporal single-cell tracing of a population of melanoma cells with a mesenchymal-like state revealed that these cells do not contribute to primary tumour growth but, instead, constitute a pool of metastatic initiating cells that switch cell identity while disseminating to secondary organs. Our data provide a spatially and temporally resolved map of the diversity and trajectories of melanoma cell states and suggest that the ability to support growth and metastasis are limited to distinct pools of cells. The observation that these phenotypic competencies can be dynamically acquired after exposure to specific niche signals warrant the development of therapeutic strategies that interfere with the cancer cell reprogramming activity of such microenvironmental cues.

Comparative single-cell RNA-sequencing profiling of BMP4-treated primary glioma cultures reveals therapeutic markers.

BACKGROUND: Glioblastoma (GBM) is the most aggressive primary brain tumor. Its cellular composition is very heterogeneous, with cells exhibiting stem-cell characteristics (GSCs) that co-determine therapy resistance and tumor recurrence. Bone Morphogenetic Protein (BMP)-4 promotes astroglial and suppresses oligodendrocyte differentiation in GSCs, processes associated with superior patient prognosis. We characterized variability in cell viability of patient-derived GBM cultures in response to BMP4 and, based on single-cell transcriptome profiling, propose predictive positive and early-response markers for sensitivity to BMP4. METHODS: Cell viability was assessed in 17 BMP4-treated patient-derived GBM cultures. In two cultures, one highly-sensitive to BMP4 (high therapeutic efficacy) and one with low-sensitivity, response to treatment with BMP4 was characterized. We applied single-cell RNA-sequencing, analyzed the relative abundance of cell clusters, searched for and identified the aforementioned two marker types, and validated these results in all 17 cultures. RESULTS: High variation in cell viability was observed after treatment with BMP4. In three cultures with highest sensitivity for BMP4, a substantial new cell subpopulation formed. These cells displayed decreased cell proliferation and increased apoptosis. Neuronal differentiation was reduced most in cultures with little sensitivity for BMP4. OLIG1/2 levels were found predictive for high sensitivity to BMP4. Activation of ribosomal translation (RPL27A, RPS27) was up-regulated within one day in cultures that were very sensitive to BMP4. CONCLUSION: The changes in composition of patient-derived GBM cultures obtained after treatment with BMP4 correlate with treatment efficacy. OLIG1/2 expression can predict this efficacy, and upregulation of RPL27A and RPS27 are useful early-response markers.

Decoding gene regulation in the fly brain.

The Drosophila brain is a frequently used model in neuroscience. Single-cell transcriptome analysis1-6, three-dimensional morphological classification7 and electron microscopy mapping of the connectome8,9 have revealed an immense diversity of neuronal and glial cell types that underlie an array of functional and behavioural traits in the fly. The identities of these cell types are controlled by gene regulatory networks (GRNs), involving combinations of transcription factors that bind to genomic enhancers to regulate their target genes. Here, to characterize GRNs at the cell-type level in the fly brain, we profiled the chromatin accessibility of 240,919 single cells spanning 9 developmental timepoints and integrated these data with single-cell transcriptomes. We identify more than 95,000 regulatory regions that are used in different neuronal cell types, of which 70,000 are linked to developmental trajectories involving neurogenesis, reprogramming and maturation. For 40 cell types, uniquely accessible regions were associated with their expressed transcription factors and downstream target genes through a combination of motif discovery, network inference and deep learning, creating enhancer GRNs. The enhancer architectures revealed by DeepFlyBrain lead to a better understanding of neuronal regulatory diversity and can be used to design genetic driver lines for cell types at specific timepoints, facilitating their characterization and manipulation.

Co-invalidation of Prnp and Sprn in FVB/N mice affects reproductive performances and highlight complex biological relationship between PrP and Shadoo.

Shadoo and PrP belongs to the same protein family, whose biological function remains poorly understood. Previous experiments reported potential functional redundancies or antagonisms between these two proteins, depending on the tissue analysed. While knockdown experiments suggested the requirement of Shadoo in the absence of PrP during early mouse embryogenesis, knockout ones, on the contrary, highlighted little impact, if any, of the double-knockout of these two loci. In the present study, we reinvestigated the phenotype associated with the concomitant knockout of these two genes using newly produced FVB/N Sprn knockout mice. In this genetic background, the combined two genes' knockout induces intra-uterine growth retardations, likely resulting from placental failures highlighted by transcriptomic analyses that revealed potential redundant or antagonist roles of these two proteins in different developmental-related pathways. It also induced an increased perinatal-lethality and ascertained the role of these two loci in the lactation process.

Drosophila Voltage-Gated Sodium Channels Are Only Expressed in Active Neurons and Are Localized to Distal Axonal Initial Segment-like Domains.

In multipolar vertebrate neurons, action potentials (APs) initiate close to the soma, at the axonal initial segment. Invertebrate neurons are typically unipolar with dendrites integrating directly into the axon. Where APs are initiated in the axons of invertebrate neurons is unclear. Voltage-gated sodium (NaV) channels are a functional hallmark of the axonal initial segment in vertebrates. We used an intronic Minos-Mediated Integration Cassette to determine the endogenous gene expression and subcellular localization of the sole NaV channel in both male and female Drosophila, para Despite being the only NaV channel in the fly, we show that only 23 ± 1% of neurons in the embryonic and larval CNS express para, while in the adult CNS para is broadly expressed. We generated a single-cell transcriptomic atlas of the whole third instar larval brain to identify para expressing neurons and show that it positively correlates with markers of differentiated, actively firing neurons. Therefore, only 23 ± 1% of larval neurons may be capable of firing NaV-dependent APs. We then show that Para is enriched in an axonal segment, distal to the site of dendritic integration into the axon, which we named the distal axonal segment (DAS). The DAS is present in multiple neuron classes in both the third instar larval and adult CNS. Whole cell patch clamp electrophysiological recordings of adult CNS fly neurons are consistent with the interpretation that Nav-dependent APs originate in the DAS. Identification of the distal NaV localization in fly neurons will enable more accurate interpretation of electrophysiological recordings in invertebrates.SIGNIFICANCE STATEMENT The site of action potential (AP) initiation in invertebrates is unknown. We tagged the sole voltage-gated sodium (NaV) channel in the fly, para, and identified that Para is enriched at a distal axonal segment. The distal axonal segment is located distal to where dendrites impinge on axons and is the likely site of AP initiation. Understanding where APs are initiated improves our ability to model neuronal activity and our interpretation of electrophysiological data. Additionally, para is only expressed in 23 ± 1% of third instar larval neurons but is broadly expressed in adults. Single-cell RNA sequencing of the third instar larval brain shows that para expression correlates with the expression of active, differentiated neuronal markers. Therefore, only 23 ± 1% of third instar larval neurons may be able to actively fire NaV-dependent APs.

Robust gene expression programs underlie recurrent cell states and phenotype switching in melanoma.

Melanoma cells can switch between a melanocytic and a mesenchymal-like state. Scattered evidence indicates that additional intermediate state(s) may exist. Here, to search for such states and decipher their underlying gene regulatory network (GRN), we studied 10 melanoma cultures using single-cell RNA sequencing (RNA-seq) as well as 26 additional cultures using bulk RNA-seq. Although each culture exhibited a unique transcriptome, we identified shared GRNs that underlie the extreme melanocytic and mesenchymal states and the intermediate state. This intermediate state is corroborated by a distinct chromatin landscape and is governed by the transcription factors SOX6, NFATC2, EGR3, ELF1 and ETV4. Single-cell migration assays confirmed the intermediate migratory phenotype of this state. Using time-series sampling of single cells after knockdown of SOX10, we unravelled the sequential and recurrent arrangement of GRNs during phenotype switching. Taken together, these analyses indicate that an intermediate state exists and is driven by a distinct and stable 'mixed' GRN rather than being a symbiotic heterogeneous mix of cells.

Cross-species analysis of enhancer logic using deep learning.

Deciphering the genomic regulatory code of enhancers is a key challenge in biology because this code underlies cellular identity. A better understanding of how enhancers work will improve the interpretation of noncoding genome variation and empower the generation of cell type-specific drivers for gene therapy. Here, we explore the combination of deep learning and cross-species chromatin accessibility profiling to build explainable enhancer models. We apply this strategy to decipher the enhancer code in melanoma, a relevant case study owing to the presence of distinct melanoma cell states. We trained and validated a deep learning model, called DeepMEL, using chromatin accessibility data of 26 melanoma samples across six different species. We show the accuracy of DeepMEL predictions on the CAGI5 challenge, where it significantly outperforms existing models on the melanoma enhancer of IRF4 Next, we exploit DeepMEL to analyze enhancer architectures and identify accurate transcription factor binding sites for the core regulatory complexes in the two different melanoma states, with distinct roles for each transcription factor, in terms of nucleosome displacement or enhancer activation. Finally, DeepMEL identifies orthologous enhancers across distantly related species, where sequence alignment fails, and the model highlights specific nucleotide substitutions that underlie enhancer turnover. DeepMEL can be used from the Kipoi database to predict and optimize candidate enhancers and to prioritize enhancer mutations. In addition, our computational strategy can be applied to other cancer or normal cell types.

Identification of genomic enhancers through spatial integration of single-cell transcriptomics and epigenomics.

Single-cell technologies allow measuring chromatin accessibility and gene expression in each cell, but jointly utilizing both layers to map bona fide gene regulatory networks and enhancers remains challenging. Here, we generate independent single-cell RNA-seq and single-cell ATAC-seq atlases of the Drosophila eye-antennal disc and spatially integrate the data into a virtual latent space that mimics the organization of the 2D tissue using ScoMAP (Single-Cell Omics Mapping into spatial Axes using Pseudotime ordering). To validate spatially predicted enhancers, we use a large collection of enhancer-reporter lines and identify ~ 85% of enhancers in which chromatin accessibility and enhancer activity are coupled. Next, we infer enhancer-to-gene relationships in the virtual space, finding that genes are mostly regulated by multiple, often redundant, enhancers. Exploiting cell type-specific enhancers, we deconvolute cell type-specific effects of bulk-derived chromatin accessibility QTLs. Finally, we discover that Prospero drives neuronal differentiation through the binding of a GGG motif. In summary, we provide a comprehensive spatial characterization of gene regulation in a 2D tissue.

The Prion-like protein Shadoo is involved in mouse embryonic and mammary development and differentiation.

Shadoo belongs to the prion protein family, an evolutionary conserved and extensively studied family due to the implication of PrP in Transmissible Spongiform Encephalopathies. However, the biological function of these genes remains poorly understood. While Sprn-knockdown experiments suggested an involvement of Shadoo during mouse embryonic development, Sprn-knockout experiments in 129Pas/C57BL/6J or 129Pas/FVB/NCr mice did not confirm it. In the present study, we analyzed the impact of Sprn gene invalidation in a pure FVB/NJ genetic background, using a zinc finger nuclease approach. The in-depth analysis of the derived knockout transgenic mice revealed a significant increase in embryonic lethality at early post-implantation stages, a growth retardation of young Sprn-knockout pups fed by wild type mice and a lactation defect of Sprn-knockout females. Histological and transcriptional analyses of knockout E7.5 embryos, E14.5 placentas and G7.5 mammary glands revealed specific roles of the Shadoo protein in mouse early embryogenesis, tissue development and differentiation with a potential antagonist action between PrP and Shadoo. This study thus highlights the entanglement between the proteins of the prion family.

A Single-Cell Transcriptome Atlas of the Aging Drosophila Brain.

The diversity of cell types and regulatory states in the brain, and how these change during aging, remains largely unknown. We present a single-cell transcriptome atlas of the entire adult Drosophila melanogaster brain sampled across its lifespan. Cell clustering identified 87 initial cell clusters that are further subclustered and validated by targeted cell-sorting. Our data show high granularity and identify a wide range of cell types. Gene network analyses using SCENIC revealed regulatory heterogeneity linked to energy consumption. During aging, RNA content declines exponentially without affecting neuronal identity in old brains. This single-cell brain atlas covers nearly all cells in the normal brain and provides the tools to study cellular diversity alongside other Drosophila and mammalian single-cell datasets in our unique single-cell analysis platform: SCope (http://scope.aertslab.org). These results, together with SCope, allow comprehensive exploration of all transcriptional states of an entire aging brain.

Mutated but Not Deleted Ovine PrP(C) N-Terminal Polybasic Region Strongly Interferes with Prion Propagation in Transgenic Mice.

UNLABELLED: Mammalian prions are proteinaceous infectious agents composed of misfolded assemblies of the host-encoded, cellular prion protein (PrP). Physiologically, the N-terminal polybasic region of residues 23 to 31 of PrP has been shown to be involved in its endocytic trafficking and interactions with glycosaminoglycans or putative ectodomains of membrane-associated proteins. Several recent reports also describe this PrP region as important for the toxicity of mutant prion proteins and the efficiency of prion propagation, both in vitro and in vivo. The question remains as to whether the latter observations made with mouse PrP and mouse prions would be relevant to other PrP species/prion strain combinations given the dramatic impact on prion susceptibility of minimal amino acid substitutions and structural variations in PrP. Here, we report that transgenic mouse lines expressing ovine PrP with a deletion of residues 23 to 26 (KKRP) or mutated in this N-terminal region (KQHPH instead of KKRPK) exhibited a variable, strain-dependent susceptibility to prion infection with regard to the proportion of affected mice and disease tempo relative to findings in their wild-type counterparts. Deletion has no major effect on 127S scrapie prion pathogenesis, whereas mutation increased by almost 3-fold the survival time of the mice. Deletion marginally affected the incubation time of scrapie LA19K and ovine bovine spongiform encephalopathy (BSE) prions, whereas mutation caused apparent resistance to disease. IMPORTANCE: Recent reports suggested that the N-terminal polybasic region of the prion protein could be a therapeutic target to prevent prion propagation or toxic signaling associated with more common neurodegenerative diseases such as Alzheimer's disease. Mutating or deleting this region in ovine PrP completes the data previously obtained with the mouse protein by identifying the key amino acid residues involved.

PlcRa, a new quorum-sensing regulator from Bacillus cereus, plays a role in oxidative stress responses and cysteine metabolism in stationary phase.

We characterized a new quorum-sensing regulator, PlcRa, which is present in various members of the B. cereus group and identified a signaling heptapeptide for PlcRa activity: PapRa(7). We demonstrated that PlcRa is a 3D structural paralog of PlcR using sequence analysis and homology modeling. A comparison of the transcriptomes at the onset of stationary phase of a ΔplcRa mutant and the wild-type B. cereus ATCC 14579 strain showed that 68 genes were upregulated and 49 genes were downregulated in the ΔplcRa mutant strain (>3-fold change). Genes involved in the cysteine metabolism (putative CymR regulon) were downregulated in the ΔplcRa mutant strain. We focused on the gene with the largest difference in expression level between the two conditions, which encoded -AbrB2- a new regulator of the AbrB family. We demonstrated that purified PlcRa bound specifically to the abrB2 promoter in the presence of synthetic PapRa(7), in an electrophoretic mobility shift assay. We further showed that the AbrB2 regulator controlled the expression of the yrrT operon involved in methionine to cysteine conversion. We found that the ΔplcRa mutant strain was more sensitive to hydrogen peroxide- and disulfide-induced stresses than the wild type. When cystine was added to the culture of the ΔplcRa mutant, challenged with hydrogen peroxide, growth inhibition was abolished. In conclusion, we identified a new RNPP transcriptional regulator in B. cereus that activated the oxidative stress response and cysteine metabolism in transition state cells.

Efficient gene expression profiling of laser-microdissected melanoma metastases.

Comparing the transcriptomes of primary and metastatic tumour tissues is a useful strategy for studying tumour progression. One factor limiting the interpretation of tissue-based transcriptomic data is the lack of cell-type purity. Laser capture microdissection (LCM) has been shown to be useful for overcoming this limitation. We established an efficient protocol for gene expression profiling of LCM and matched metastatic melanomas using a transgenic mouse model. This optimized workflow combines microsurgical recovery of mouse lungs, appropriate tissue freezing, laser microdissection of homogeneous tumour cell populations from cryosections, isolation of high-quality RNA and gene expression analysis. The RNA isolated from laser-microdissected material was not contaminated by stroma cells, was of excellent quality, and the synthesis of cDNAs was homogeneous and highly reproducible. Subsequent custom-based Taqman-low-density-array (TLDA)-based gene expression profiling identified stronger expression of five genes (M-MITF, TYR, STAT3, CCND1 and PAX3) in primary than metastatic melanoma. We detected only minor transcriptomic differences between primary and metastatic melanoma tissue. This optimized workflow could be very valuable for various studies requiring cell type-specific transcriptomic analysis.

Prion protein and Shadoo are involved in overlapping embryonic pathways and trophoblastic development.

The potential requirement of either the Prion or Shadoo protein for early mouse embryogenesis was recently suggested. However, the current data did not allow to precise the developmental process that was affected in the absence of both proteins and that led to the observed early lethal phenotype. In the present study, using various Prnp transgenic mouse lines and lentiviral vectors expressing shRNAs that target the Shadoo-encoding mRNA, we further demonstrate the specific requirement of at least one of these two PrP-related proteins at early developmental stages. Histological analysis reveals developmental defect of the ectoplacental cone and important hemorrhage surrounding the Prnp-knockout-Sprn-knockdown E7.5 embryos. By restricting the RNA interference to the trophoblastic cell lineages, the observed lethal phenotype could be attributed to the sole role of these proteins in this trophectoderm-derived compartment. RNAseq analysis performed on early embryos of various Prnp and Sprn genotypes indicated that the simultaneous down-regulation of these two proteins affects cell-adhesion and inflammatory pathways as well as the expression of ectoplacental-specific genes. Overall, our data provide biological clues in favor of a crucial and complementary embryonic role of the prion protein family in Eutherians and emphasizes the need to further evaluate its implication in normal and pathological human placenta biology.

Toxic lipids stored by Kupffer cells correlates with their pro-inflammatory phenotype at an early stage of steatohepatitis.

BACKGROUND & AIMS: Non-alcoholic steatohepatitis (NASH) is characterized by steatosis associated with liver inflammation. Steatosis causes recruitment of lymphocytes into the liver and this is worsened by lipopolysaccharides (LPS). As macrophages may be involved in the lymphocyte homing, we studied the role of lipids in determining the phenotype of Kupffer cells (KCs) at the stage of steatosis. METHODS: Steatosis was induced in mice by a high fat diet. The turnover and the recruitment of KCs were analyzed in vivo by flow cytometry. KCs phenotype was assessed by optical and electron microscopy, cell culture and lymphocyte recruitment by in vitro chemotaxis. Lipidomic analysis was carried out by mass-spectrometry and gene expression analysis by TaqMan low density array. RESULTS: Although the number of KCs was not modified in steatotic livers compared to normal livers, their phenotypes were different. Electron microscopy demonstrated that the KCs from fatty livers were enlarged and loaded with lipid droplets. Lipid synthesis and trafficking were dysregulated in fat-laden KCs and toxic lipids accumulated. Fat-laden KCs recruited more CD4+ T and B lymphocytes in response to LPS stimulation than did control KCs and produced high levels of pro-inflammatory cytokines/chemokines, which could be reversed by inhibition of lipogenesis. CONCLUSIONS: Lipid accumulation in fat-laden KCs is due to a dysregulation of lipid metabolism and trafficking. Fat-laden KCs are "primed" to recruit lymphocytes and exhibit a pro-inflammatory phenotype, which is reversible with inhibition of lipogenesis.

Maintaining RNA integrity in a homogeneous population of mammary epithelial cells isolated by Laser Capture Microdissection.

BACKGROUND: Laser-capture microdissection (LCM) that enables the isolation of specific cell populations from complex tissues under morphological control is increasingly used for subsequent gene expression studies in cell biology by methods such as real-time quantitative PCR (qPCR), microarrays and most recently by RNA-sequencing. Challenges are i) to select precisely and efficiently cells of interest and ii) to maintain RNA integrity. The mammary gland which is a complex and heterogeneous tissue, consists of multiple cell types, changing in relative proportion during its development and thus hampering gene expression profiling comparison on whole tissue between physiological stages. During lactation, mammary epithelial cells (MEC) are predominant. However several other cell types, including myoepithelial (MMC) and immune cells are present, making it difficult to precisely determine the specificity of gene expression to the cell type of origin. In this work, an optimized reliable procedure for producing RNA from alveolar epithelial cells isolated from frozen histological sections of lactating goat, sheep and cow mammary glands using an infrared-laser based Arcturus Veritas LCM (Applied Biosystems®) system has been developed. The following steps of the microdissection workflow: cryosectioning, staining, dehydration and harvesting of microdissected cells have been carefully considered and designed to ensure cell capture efficiency without compromising RNA integrity. RESULTS: The best results were obtained when staining 8 µm-thick sections with Cresyl violet® (Ambion, Applied Biosystems®) and capturing microdissected cells during less than 2 hours before RNA extraction. In addition, particular attention was paid to animal preparation before biopsies or slaughtering (milking) and freezing of tissue blocks which were embedded in a cryoprotective compound before being immersed in isopentane. The amount of RNA thus obtained from ca.150 to 250 acini (300,000 to 600,000 µm2) ranges between 5 to 10 ng. RNA integrity number (RIN) was ca. 8.0 and selectivity of this LCM protocol was demonstrated through qPCR analyses for several alveolar cell specific genes, including LALBA (α-lactalbumin) and CSN1S2 (αs2-casein), as well as Krt14 (cytokeratin 14), CD3e and CD68 which are specific markers of MMC, lymphocytes and macrophages, respectively. CONCLUSIONS: RNAs isolated from MEC in this manner were of very good quality for subsequent linear amplification, thus making it possible to establish a referential gene expression profile of the healthy MEC, a useful platform for tumor biomarker discovery.

Enterococcal leucine-rich repeat-containing protein involved in virulence and host inflammatory response.

Enterococcus faecalis is an important nosocomial pathogen associated with high morbidity and mortality for patients who are immunocompromised or who have severe underlying diseases. The E. faecalis genome encodes numerous surface-exposed proteins that may be involved in virulence. This work describes the characterization of the first internalin-like protein in E. faecalis, ElrA, belonging to the recently identified WxL family of surface proteins. ElrA contains an N-terminal signal peptide for export, a leucine-rich repeat domain that may interact with host cells, and a C-terminal WxL domain that interacts with the peptidoglycan. Disruption of the elrA gene significantly attenuates bacterial virulence in a mouse peritonitis model. The elrA deletion mutant also displays a defect in infection of host macrophages and a decreased interleukin-6 response in vivo. Finally, elrA expression is induced in vivo. Altogether, these results demonstrate a role for ElrA in the E. faecalis infectious process in vivo and suggest that this surface protein may contribute to E. faecalis virulence by stimulating the host inflammatory response.